Wave-packet dynamics in slowly perturbed crystals: Gradient corrections 

and Berry-phase effects 
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We present a unified theory for wave-packet dynamics of electrons in crystals subject to pertur- 
bations varying slowly in space and time. We derive the wave-packet energy up to the first order 
gradient correction and obtain all kinds of Berry-phase terms for the semiclassical dynamics and 
the quantization rule. For electromagnetic perturbations, we recover the orbital magnetization en- 
ergy and the anomalous velocity purely within a single-band picture without invoking inter-band 
couplings. For deformations in crystals, besides a deformation potential, we obtain a Berry-phase 
term in the Lagrangian due to lattice tracking, which gives rise to new terms in the expressions 
for the wave-packet velocity and the semiclassical force. For multiple-valued displacement fields 
surrounding dislocations, this term manifests as a Berry phase, which we show to be proportional 
to the Burgers vector around each dislocation. 

PACS numbers: 72.15.-v, 72.10.Di, 72.10.Fk 



I. INTRODUCTION 

Our understanding of electronic properties of crys- 
talline solids is pjijimarily based on the Bloch theory for 
periodic systems^ It has been of great interest to extend 
this theory to situations where crystals are perturbed in 
various ways. So far the most useful description has been 
the semiclassical theory for electron dynamics within a 
band supplemented by the semiclassical quantization rule 
or the Boltzmann transport equations. For example, the 
equations of moticui of Bloch electrons in electromagnetic 
fields are given byo 

1 d£oAk) 



dk ' 
?ik = — eE — ex X B, 



(1.1) 



where fo,n(k) is the energy of the nth band of an un- 
perturbed crystal. These equations have played a funda- 
mental role in the physics of metals and semiconductors. 

The derivation of Eq. (1.1) dates back to Bloch, Peierls, 
Jones and Zener in the early 1930'sB By assuming that 
the transition probabilities to o ther bands are negligible, 
they showed that the Eqs. (1.1) describe the motion of a 
narrow wave packet obtained by superposing the Bloch 
states of a band. Various extensions of the theory have 
been made to deal with perturbations of more general 
nature and to obtain corrections to Eqs. (1.1) in high 
fields. |-| 

PeierlsQ pioneered the effort of constructing an effective 
one-band Hamiltonian to describe the quantum dynamics 
of a Bloch electron. By using the tight-binding model, 
he was able to show that the effective Hamiltonian in 
the presence of a magnetic field may be obtained by re- 
placing the crystal momentum hk. by the gauge invariant 
momentum operator [—iK\7 + eA(x)] in the unperturbed 
band energy: 



^^off — £o,n 



-iV + ^A(x) 
n 



(1.2) 



which later came to be known j-as the Peierls substitu- 
tion. Two decades later, SlateiEl and LuttingerEl gave a 
more rigorous derivation of the effective Hamiltonian for 
electromagnetic perturbations, by expanding the wave 
function in the basis of Wannier functions: 



^[/(x,^ =^/z(t)M^(x 



(1.3) 



where {R;} are the lattice positions. They showed that 
the envelope function /(x, i), defined by /(R;,i) ~ fi{t) 
and a smooth interpolation between the atomic positions, 
satisfies the effective Schrodinger equation 

in-^J = [£o,n + Ja(x)] - e(/)(x)} /, (1.4) 

where 0(x) is a slowly varying scalar potential. The equa- 
tions of motion (1.1) then follow from Eq. (1.4) and the 



correspondence principle. 

Further development of the theory was made by taking 
into account the effects of inter-band coupling. AdamsQ 
extended the works of Slater and Luttinger to many band 
operator formalism. Karplus, Luttinger, and Kohn de- 
rived a correction to the velocity, known as the anoma- 
lous velocity, and predicted a spontaneous Hall effect in. 
ferromagnetic materials.cl Later, Adams and BlounmH 
showed that this term arises from the noncommutabil- 
ity between the Cartesian components of th e ipt ra-band 
position operator. Recently, Chang and Niullililj related 
the anomalous velocity correction to the Berry phas^as-, 
sociatcd with the electron motion in an energy band.t3~£3 
Corrections to the effective Hamiltonian as an asymptotic 
series in the field strength were obtained by eliminating 



the inter-band matrix elements with unitary transfornxELi 
tions by Kohp- Blount, and Roth in the early 1960's.ll3 
Later, BrownliZI extended the Wannier function method 
to crystals under high magnetic fields using magnetic 
translation symmetry. A decade later, a variant of the 
Wannier function method that treats position and mo- 
mentum in a symmetriG-way, known as fcg-representation, 
was developed by Zak.H Recenldy, Rammal and Bellis- 
sard used an alaebraic approachllj and Wilkinson an op- 
erator approachcJ to derive the first order field correction 
for the special case of the Harper's equation.Ljl£3 

Another approach to this problem involves the use of 
the WKB expansion to derive a Hamilton- Jacobi equa- 
tion, and then making the correspondence from the clas- 
sical variables to the quantum operators. This method 
was applied by many researchers to understand the na- 
ture of the spectrum and the wave function of electrons 
described by the Ilarper's equation, particularly in the 
semiclassical limit£3 A more general treatment of the 
problem is based on a two-scale expansion in which the 
electron coordinate and the slowly varyingjjzjsctor poten- 
tial are regarded as independent variables .Ej 

In this work we come back to the original considera- 
tion of a wave packet in ^a, band and use a time depen- 
dent variational principlaHJ'E^ to derive the wave-packet 
energy up to first order in the gradient of the pertur- 
bations and Berry-phase corrections to the semiclassical 
dynamics and the quantization rule. We are able to ob- 
tain the magnetization energy and the anomalous veloc- 
ity entirely from the single-band point of view. Also, our 
method can be directly extended to the case of slowly 
perturbed magnetic bands, where methods based on the 
usual Wannier function approach break down because 
of nonexistence oL localized Wannier functions for such 
bands in generalHj nn 

This program was started with Chang and NiuU'll3 for 
the special case of a two-dimensional periodic system in a 
strong magnetic field. Here we establish a unified frame- 
work for slowly perturbed crystals whose Hamiltonian 
can be expressed in the form 

i/[x,p;/3i(x,i),...,/3,(x,i)], (1.5) 

where {/3i(x, f)} are the modulation functions character- 
izing the perturbations. They may represent either defor- 
mation strain fields, gauge potentials of electromagnetic 
fields, or slowly varying impurity potentials. They also 
appear in modeJj)otentials for modulated and inccmmen- 
surate crystals,tZl and for graded semiconductors.E3 

We shall illustrate our formalism with two special cases 
of perturbations: electromagnetic fields and deforma- 
tions in crystals. In the first case, in addition to the 
corrections of orbital magnetization energy and anoma- 
lous velocity mentioned above, we discuss the Peierls sub- 
stitution. Berry-phase-modified Landau levels, and Zak- 
phase-modified Wannier-Stark levels. For deformational 
perturbations, we show that the deformation correction 
to the wave-packet energy can be obtained from the dif- 
ferential shift in the band energy under uniform strain. 



Then we obtain for the Lagrangian a Berry-phase term 
due to lattice tracking, which gives rise to new terms in 
the expressions for the wave-packet velocity and the semi- 
classical force; for multiple- valued displacement fields in 
the presence of dislocations, this term manifests as a 
Berry phase, which we show to be proportional to the 
Burgers vector around each dislocation, and thus in a 
sense keeps track of the lattice position. We also dis- 
cuss the consequences of the Berry-phase term on elec- 
tron transport, and the Aharonov-Bohm type effects in 
dislocated crystals. 

The paper is organized as follows. We present our for- 
malism in Sec. ^ treat ele ctro magnetic and deforma- 
tional perturbations in Sees. |lll| and |^ respectively, and 
conclude with a summary in Sec. M. 



II. FORMALISM 

We shall begin by constructing a basis local to the 
wave packet and describe the wave packet in detail. Then 
we derive the Lagrangian, the semiclassical equations of 
motion, and Berry-phase correction to the semiclassical 
quantization, and discuss some aspects of formal quanti- 
zation through the Hamiltonian formalism. 



A. The local basis and the wave packet 

Consider a wave packet centered at Xc at a given time, 
with its spread small compared to the length scale of the 
perturbations. Then the approximate Hamiltonian that 
the wave packet "feels" may be obtained by linearizing 
the perturbations about the wave packet center as 

fj TT \ sum of terms oc 1 in ^\ 

^ * + \ (i - X,) . grad,^ A(x„ t) ] ' ^^-^^ 

where He = i?[x, p; {/9i(xc, i)}] will be called the local 
Hamiltonian. The terms within the braces are small in 
the neighborhood of the wave packet and may be treated 
perturbatively. The local Hamiltonian He has the re- 
quired periodicity of the unperturbed crystal, and has 
an energy spectrum of bands (Bloch bands) with Bloch 
wave eigenstates satisfying 

-ffc(xc,i)|V'q(x„i)) =f,(x„q,i)|V'q(x„i)), (2.2) 

where q is the Bloch wave vector and fc(xc,q, t) is the 
band energy. Since we will be concerned with only a sin- 
gle band, we have omitted the band index for simplicity 
of notation. We note that both the wave-packet center 
Xc and time t enters in the Bloch states and the band 
energy parametrically. We shall see that the dependency 
on the center position of the wave packet will manifest 
as new types of Berry-phase terms in the equations of 
motion. 



These eigenstates form a convenient basis to expand 
the wave packet. Specifically, we write 

1^-)= y"Aa(q,t)|^q(x„i)), (2.3) 

where a(q, t) is the amplitude with the normalization 

d\|a(q,t)|2-(*|*) = l. (2.4) 



(2.5) 



Here we have used the convention that 
(V'q|^q') = <5(q-q'). 



It is assumed that the distribution |a(q, t)p is narrow 
compared to the size of the Brillouin zone and has the 
mean wave vector 



Aq|a(q, t)\' 



(2.6) 



To be consistent, the wave packet must yield the pre- 
assigned center position: 



(2.7) 



This condition can be expressed in terms of other wave- 
packet parameters as follows. Writing the amplitude in 
the form a(q, t) — |a(q, t)| exp[— i7(q, i)], and using the 
matrix elements of position operator x between the Bloch 
states [Eq. (Al) of Appendix we find that 



q \a[ 



9q 



. du 
9q 



(2., 



where |u) = |M(xc,q, t)} 



— „-iqx 



|'!/'q(xc, i)) is the peri- 



odic part of the Bloch wave, and inner products involving 
the periodic part \u) mean an integration over the unit- 
cell volume Vc with a factor of {2n)'^ /vc^ which implies 
the normalization {u\u) — 1. According to our assump- 
tion of a narrow wave packet in the q-space, Eq. (2.7) 
becomes 



97c 
9qc 



du 



(2.9) 



where \u) now stands for |u(xc, qc, ^)), and 7c for 7(qc, t). 

In writing down the expansion (2.3), we have assumed 
that the wave packet that is initially in a band always lies 
in the same band. This is justified if the band is sepa- 
rated from other bands by finite gaps, and if the time and 
the length scales of the pcrturbatio]i*,are long compared 
to those associated with these gaps.c3 



B. The lagrangian and dynamics 

The dynamics of the mean position Xc and the crys- 
tal momentum ?iqc can in principle be derived from the 
Schrodinger equation for the wave packet. It is more 



convemep 
principle 



:,obtained using a time dependent variational 
with the Lagrangian given by 



L={^i--H\^) 



(2.10) 



where and hereafter we use the convention fi = 1. We 
use d/dt to mean the derivative with respect to the time 
dependence of the wave function explicitly or implicitly 
through Xc and qc. The partial derivative d/dt, is re- 
served for those with Xc and qc held fixed. 

Under the previously discussed conditions on the 
widths of the wave packet, we can evaluate the La- 
grangian as a function of Xc and qc, their time deriva- 
tives, and the time t: 



L « L(Xc,Xc,qc,qc,i)- 



(2.11) 



The terms involving higher moments of the wave packet, 
which specify its width and shape, are of higher or- 
der in thfi-ffradient of the perturbations and hence are 
neglected. E3 



Accordingly, we obtain for the first term in Eq. (2.10) 



.d^ 

'~dt 



dt 



.du 



The first term comes from the explicit time dependence 
of 7c. The contribution fro m |Q (q, t)\ is zero because of 
the normalization condition (2.4) on the amplitude. The 
second and third terms come about because of the de- 
pendence of the basis functions on time explicitly and 
implicitly through Xc respectively. Using the relation 



dt 



dt 



9qc 



(2.13) 



and Eq. (2J), Eq. ( 2.12 ) can be recast into the form 



~dt 



d'fc 
dt 



. du 

'9qc 

du\ 
^^), (2.14) 



where 7c appears only in a total time derivative. 

The expectation value of the Hamiltonian, which gives 
the wave-packet energy f , may be evaluated up to first 
order in the perturbation gradients using the linearized 
Hamiltonian (2.1): 



£ = ~ (*|i?c|*) + (*|Ai7|*}, (2.15) 

where the gradient correction Ai/ may be written aj^^ 



-4 



(X - Xc 



dHr dHr 



dy 



dy 



(X - Xc) 



(2.16) 



The expectation value of the local Hamiltonian is just 
the band energy at the mean wave vector, 



(^-IFcl*) =£c(xc,qc,i), 



(2.17) 



C. Formal and semiclassical quantization 



while the gradient correction requires some calculations 
[Appendix Q , but the result has the simple form 



-Im 



du 



{8c — He 



du 



(2.18) 



where "•" denotes scalar product between the vectors 
formed by gradients with respect to Xc and q. 
The Lagrangian thus takes the form 



L= -£ 



+ Xc 



qc ■ Xc + qc • \ u 

. du 
9xc 



, du 

du\ 

'ml 



(2.19) 



where we have neglected a term of total time deriva- 
tive d(7 — Xc • (\c)/dt in the Lagrangian, as it does not 
affect the equations of motion and the quantization 
rule. The last three terms may be grouped into a sin- 
gle term {u\idu/dt), which turns out to be the net 
rate of change of Berry phase for wave-packet motion 
within the band. We note that under the transforma- 
tions \u) — > exp[i(/?(xc, q, t)] \u) or q ^ q + K, K being a 
reciprocal lattice vector, the Lagrangian remains invari- 
ant up to a total time derivative of some function of Xc, 
qc, and t. The former corresponds to gauge invariance 
while the latter to periodicity in the reciprocal space. 

From the Lagrangian ( 2.19| ) we obtain the following 
equations of semiclassical motion: 

d£ 

dqc 
d£ 

qc = h (fixx • Xc + rixq ■ qc) - rJtx- (2.20) 



The components of the tensor fJqq are defined by 



Go ^ 
Xc — (^^qx ■ Xc -t- ^qq ' qc J ~t" ^tq? 



(f^qq)Q/3 = ^ 



9a 9/3 



du 



dqc 



du 



dq, 



cf3 



di 



dq 



du 



dqc 



(2.21) 



and those of the vector fttx by 

du 



du 
'dt 



dxc 



du 



dxc 



du 
'dt 



(2.22) 



where a and /3 are Cartesian indices. The other ten- 
sors rixx, f^xq, and fiqx and the vector n^q are de- 
fined similarhr,. These quantities are known as Berry 
curvatures .E3o We note that these equations involve 
Berry curvatures between every pair of parameters and 
that they have symplectic symmetry in the absence of 
time dependence. 



Wc mentioned in the Introduction that the equations 
of motion were usually derived from the effective Hamil- 
tonian upon using the correspondence principle. Here 
we consider the reverse process to obtain the effective 
quantum Hamiltonian from the semiclassical dynamics. 
This requires a knowledge of canonical structure of the 
wave-packet dynamics. Following the standard proce- 
dure of analytical mechanics, we introduce the canonical 
momenta conjugate to the generalized coordinates 



Pi = — 
i9xc 

dqc 



u 

. du 
dqc 



du 

^ 

dxc 



(2.23) 
(2.24) 



and the semiclassical Hamiltonian Ti. by the Legendre 
transformation 



7^ = Xc • Pi + qc • P2 - 



£(Xc,qc,i) 



, du 



(2.25) 



The semiclassical Hamiltonian is independent of Pi and 
P2, because the Lagrangian is linear in the generalized 
velocities. 

Starting with Ti., regarded formally as a function of Xc, 
qc, Pi, and P2, one cannot obtain the equations of mo- 



tion ( 2.20|) from the H amilton equations. This is because 
Eqs. ( 2.23| ) and ( ^.24| ) defining Pi and P2 do not depend 
on the generalized velocities, and hence they should be 
treated as constraints between the canonical variables. 

In the simple case where the Berry-phase terms are 
zero, these constraints become Pi = qc, and P2 = 0, 
and H — £■ These suggest that we treat (xc,qc) as a 
canonical pair and forget about the other degrees of free- 
dom. By d oing so, we can indeed obtain the equations 
of motion ( ^.20 ) from the Hamilton equations. Hav- 
ing identified the canonical pair, one can proceed with 
a formal procedure of quantization ("requantization"), 
qc ~id/dxc, to obtain an effective quantum Hamilto- 
nian. A slightly more general case will be encountered 
in the case of electromagnetic perturbations in the next 
section of this article. 

Wh en B erry-p hase terms are present, constraints 
( ^.231) and ( |2.24 still imply some hidden canonical re- 
lations between Xc and qc, but these are entangled in a 
complicated manner that cannot be expressed explicitly 
in general. This clearly shows the difficulty of the Hamil- 
tonian formalism in the presence of Berry-phase terms. 
If one insists on using the Hamiltonian approach, one can 
employ the method of Lagrange multipliers, which allows 
the spurious degsees of freedom to be formally treated 
as ind ependent and obtain the equations of motion 
(^.20| ). The "requantization" procedure for this case is 
quite complicated and needs further investigation. 



The semiclassical quantization, on the other hand, is 
quite straight forward. In order that stationary states 
and energy levels can be talked about, we shall restrict 
ourselves to static perturbations. For a wave-packet mo- 
tion that is regular and is described by closed orbits in 
the phase space (xc,qc), semiclassical enecgy levels are 
obtained using the quantization procedureE3 due to Ein- 
stein, Brillouin, and Keller, 



Pi • dxc 



dqc = 27r[m+ -], (2.26) 



where C denotes an orbit of constant energy m an 
integer that labels the eigenval ue, an d ly t he nu mber of 
caustics traversed. With Eqs. (2.23) and (2.24) for Pi 
and P2, the above condition reduces to 



V T{C) 

™ + T 7^ — 

4 27r 



where 



no 



. du 



dqc 



. du 
dqc 



(2.27) 



(2.28) 



is the Berry phase acquired by the wave packet upon 
going round the closed orbit once. The Berry-phase cor- 
rection to the quantization cpadition made its first ap- 
pearance in Wilkinson's woi'kcJ on Harper's equation, in 
Kuratsuji and lida's worktj on adiabatic nuclear motion, 
and also recently in the work of Chang and Niuli3 on 
wave-packet dynamics in magnetic Bloch bands. 



III. ELECTROMAGNETIC FIELDS 



This has the form ( |l.5| ) and (3.1), with the gauge poten- 
tials playing the role of the modulation functions, and 
hence the local Hamiltonian must have the form 



He = -ffo[q + eA(xc,<)] - e(j>{xc,t). 



(3.3) 



As eA(xc, t) is only an additive constant to the crys- 
tal momentum q, the basis states have the form 
|u(xc,q, i)) = |w(k)), where k = q-|-eA(xc,i) is the 
gauge invariant or mechanical crystal momentum. In 
terms of the gauge invariant crystal momentum k, a num- 
ber of simplifications can be obtained. First, the eigenen- 
ergy can be written in the form 



fc(xc,k,t) =£:o(k) 



(Xc,t). 



(3.4) 



Second, the gradient correction ( 2.18| ) becomes the or- 
bital magnetization energy of the wave packet, 

-M-B, (3.5) 

where B = curlx^ A(xc, i) is the magnetic field, and 



M 



-e Im 



du 
dk 



(£o--ffo(k)) 



du 
dk 



(3.6) 



k=ke 



is the orbital magnetic moment of Bloch electrons. Third, 
the last three terms of the Lagrangian ( 2.19| ) simply be- 
come the single term k^ ■ {u\idu / dkc) . Finally, the La- 
grangian takes the form 



L = + e(/)(xc, t) +±c- kc 
- exc ■ A(xc,t) + kc 



du . , , 



So far our treatment of perturbations has been in gen- 
eral terms, and our results are in an abstract form. Their 
physical meaning will become clear through the consid- 
eration of two special cases in this and the next section. 
For a class of perturbations for which the Hamiltonian is 
of the special form 



Ho [x + /3i (x, t), p + /32(x, t)] + /33(x, t), 



(3.1) 



all the results can be expressed in terms of the unper- 
turbed Bloch wave basis. In this section we shall consider 
electromagnetic perturbations for which /3i(x,t) = 0. 



A. The gauge invariant crystal momentum 

Let i?o(q) denote the Hamiltonian for the bare crystal, 
with the eigenstate |u(q)) (the periodic part of the Bloch 
wave) and the band energy £'o(q) for a particular band. 
The Hamiltonian gets modified by the gauge potentials 
[A(x, t), (/)(x, t)] of an electromagnetic field to 



H = Ho[q + eA{±,t)] - ecj){±,t). 



where 8^ = £o(kc) - M ■ B. 



B. The reciprocal magnetic field and orbital 
magnetization energy 

The equations of motion can cither be derived vari- 
ation ally f rom the above Lagrangian or directly from 
Eqs. ( 2.20 ) derived for the general case in the previous 
section. They have the form 



• _ ir o 

akc 

kc = — e E — e Xc X B, 



(3.8) 



where E = — gradx^0(xc, t) — dA{xc,t)/dt is the electric 
field, and 



(3.9) 



are the components of the vector form of the antisymmet- 



(3.2) ric tensor ilkk given by (2.21). In (3.£) and henceforth. 



repeated Cartesian indices are taken to be summed. Be- 
cause rj occupies a similar position as the magnetic field 
in the equations of motion, it will be called the reciprocal 
magnetic field. 



Equations of motion (3.8) differ from Eqs. (1.1) in 



two respects. Firstly, the energy £ contains a correction 
term from the magnetic moment (3.6) of Bloch electron. 
This term has been derived earlier as a first order correc- 
tion iit-the theory of Bloch electrons subject to magnetic 
fields.Ea A similar term has also been found-in-the the- 
ory of electrons in incommensurate lattices,E3tJ and in 
the theory of wave packet dynamics in magnetic Bloch 
bands.E3 Secondly, the correction term to the velocity, 
— kc X Jl, is the anomalous velocity that was predicted 
to give rise to a spoojLaneous Hall conductivity in fer- 
romagnetic materials.aLi In the context of the quantum 
Hall effect, the integral of the Berry curvature ^k^kf, over 
the Brillouin zone was shown to be proportional to the 
Hall conductivity for a full band and to be quantized 
(Chcsn's topological invariant ).EZl Recently, Chang and 
NiutHi proved this result semiclassically. It seems more 
appropriate to call this term Hall velocity than anoma- 
lous velocity. 

should be invariant 



The semiclassical equations (3 



under time reversal, spatial inversion, or certain rotations 
if these are symmetries of the unperturbed crystal. Such 
symmetries impose severe restrictions on the behavior of 
the reciprocal magnetic field $7 and the magnetic mo- 
ment M as functions of k. Under time reversal, Xc, kc, 
and B change sign while Xc, kc, and E are invariant. If 
the bare crystal is invariant under time reversal, we must 
have f2(-k) = -Jl(k), and M(-k) = -M(k), which im- 
plies in particular that they must vanish at k = 0. Under 
spatial inversion, E, Xc, kc, and the time derivatives of 
the last two change sign while B remains unchanged. If 
the bare crystal has inversion symmetry, we must have 
n(-k) = J7(k), and M(-k) = M(k). Finally, if the sys- 
tem is invariant under certain proper rotations, the re- 
ciprocal magnetic field and the magnetic moment should 
transform like vectors under these rotations. 

For monatomic non-magnetic crystals, both time re- 
versal and spatial inversion symmetries are present, ren- 
dering Jl and M null everywhere in the Brillouin zone. 
However, it is not entirely justified to ignore these quanti- 
ties for magnetic crystals or non-magnetic crystals with- 
out inversion symmetry (such as GaAs). Investigations 
have been undertaken to see whether the presence of the 
reciprocal magnetic field and the orbital magnetization 
lead to observable effects. 



The Lagrangian (3/7) and the equatioiis,of motion (3.g) 
were derived earlier by Chang and Niuli3 for perturbed 
magnetic Bloch electrons in two dimensions in the gauge 
where 0(x, t) — 0. They used a wave packet that gauges 
away the vector potential locally so that their magnetic 
Bloch wave vector is the same as the corresponding gauge 
invariant crystal momentum q + eA(xc,t) = k here. 
Their derivation was less general in that it is only for 
two dimensions and more general in that it provides a 



description of electrons in rational magnetic fields Bg, 
for which the flux through a unit cell equals a rational 
fraction of the flux quantum {h/e). However, 

our formalism can easily be generalized to this situation, 
for a more general three dimensional case, by assuming 
a background of constant rational magnetic field Bq. All 
we need to do is to interpret q as the wave vector of 
magnetic Bloch states defined withii|i_a reduced Brillouin 
zone (the magnetic Brillouin zone),E3 and to replace B 
in Eq. (U) by B - Bo. 



C. Peierls substitution and Landau levels 



It follows from the Lagrangian (^J) that the canonical 
momenta are given by 



Pi = kc - eA(xc,t), 
du 
9k: 



(3.10) 
(3.11) 



and the Hamiltonian by H = £m0^c) — e^(xc,i). In 
the absence of the Berry-phase term ( p.ll ), we may ob- 
tain the Hamiltonian as a function of the canonical pair 
(xc,Pi) as 



H = £m [Pi + eA(xc, t)] - e0(xc, t). 



(3.12) 



The quantization of this Hamiltonian by setting Pi = 
~id/dxc amounts to the Peierls substitution. However, 
it is not clear how to deal with the case with a Berry- 
phase term using the Hamiltonian approach. 

When only a unifor m m agnetic field is present, the two 
equations of motion ( |3.8| ) can be combined into a single 
one for the k-space motion. It is evident that a k-space 
orbit must lie in a plane normal to B and must be on a 
constant energy surface of ^(k). If such an orbit is closed, 
known as a cyclotron orbit, the EBK formula yields 



lB.£kcXdkc=^ 



1 m\ 



2tt I 



where 



r(c) 



dk- 



. du 



(3.13) 



(3.14) 



is the Berry phase accumulated by the wave packet upon 
completing a circuit along the loop C, and we have 
restored the Planck constant. The left-hand side of 



Eq. (3.13) is just the k-space area enclosed by the or- 
bit C. As this phase influences energy levels, it affects 
the density of states. It is shown in Ref. 12 that F plays 
an important role in determining the spectral splitting 
pattern of magnetic bands. 



D. Zak phase and Wannier-Stark ladder 

In this subsection, we restrict our discussion to one di- 
mension for simplicity. The semiclassical motion under 
a uniform electric field is described by the Hamiltonian 
H = £o{kc) + eExc- It follows from the boundedness 
and the periodicity of the band energy that the motion 
in real space is also bounded and periodic. Such a closed 
motion in the real space is known as Bloch oscillations. 
In the reduced zone scheme, the motion is closed also in 
the phase space (xc,kc), whic h is a cylinder. Quantizing 
this motion according to Eq. (2.26) gives the condition 



Tr/a 

j dk^Xc{K) ^ 2ti [m +'^- ^ , (3.15) 



— TT /a 

where a stands for the lattice constant, and 



TT /a 

dk { u 

-TT / a 



dk 



(3.16) 



is known as the Zak phase J3 and Xc{kc) is the con- 
stant energy curve for the mth energy level defined by 
Wm = £oikc) + eExc, m being an integer. Averaging 
this expression over the orbit, we obtain from Eq. (3.15) 



£q + eEa 



r 

2^ 



(3.17) 



where is the average of the band energy over the Bril- 
louin zone, and m is any integer between — oo to oo, since 
the mean value of Xc can be anywhere on the cylinder. 
This spectrum, known aa,the Wannier-Stark ladder, was 
first derived by WannieiEj without the Berry-phase term. 
The correction was. due to Zak, who later interpreted it 
as a Berry phase. e2I 



IV. DEFORMATIONS IN CRYSTALS 

We shall now come to dcformational perturbations. It 
turns out that the model Hamiltonian for a deformed 
crystal also has the special form (B.l) with /32(x, t) = 0, 
which is i?o[x -I- /3i(x, t), p] + /33(x, t), and hence all the 



corrections are expressible in terms of the undeformed 
basis for this case too. 



A. The translated crystal basis 



where u(x) is a smooth displacement field@ satisfying 
u(R; + u;) = Ui, and Saf} = dua/dxfs is the unsym- 
metrized strain. The justification of the above Hamilto- 
nian and an explicit expression for the last term are given 
in Appendix While the last term of the Hamiltonian, 
being proportional to the strain, can be treated perturba- 
tively, the other terms are of the form of Eq. (1.5), with 
the displacement field playing the role of the modulation 
function. The local Hamiltonian is then given by 



Hc = t; ^ Vb[x - u(xc 

Zm 



(4.2) 



This is nothing but the Hamiltonian of an undeformed 
crystal shifted in position by the displacement field at 
the center of the wave packet, u(xc). The band energy 
is therefore the same as that of the undeformed crystal, 
£o(k), and the eigenstates are the translated Bloch waves 
{i/'k[x-u(xc)]}. 

Our wave packet will thus be formed out of these trans- 
lated Bloch states of the undeformed but translated crys- 
tal. This procedure is valid so long as the strain is weak, 
so that the variation of the displacement within the spa- 
tial width of the wave packet is small. When the first 
order corrections to the Hamiltonian are taken into ac- 
count, our method should give the same physical results 
(to the same order) as obtained using a strained basis. 
However, our formulation should be simpler and easier 
to interpret, because it avoids the necessity of transfor- 
mation between the lab and lattice frames of reference. 

Although the small strain regime covers the vast ma- 
jority of practical situations, it is some times necessary 
to consider the effect of large strains. In the case where a 
large uniform and static strain is superposed on top of a 
small varying strain, our formulation can still be applied; 
one only needs to interpret the basis as that of the uni- 
formly strained crystal. There can be a third possibility 
in which the strain variation is large over large distances 
but is small over the size of the wave packet. In this case, 
it is more appropriate to use a strained local basis, that 
is, a basis of a homogeneously strained crystal with the 
strain value given by the actual strain at the center of 
the wave packet. 



B. The crystal deformation potential 

The wave-packet energy is obtained by summing the 
expectation values of th e local Hamiltonian ( |4.2| ), the 
gradient correction ( 2.1(^ ), and the last term of Eq. (p|). 
Because of the functional form of the basis states, we 
have 



A deformed crystal with atomic displacements {u;} 
may be described by the Hamiltoniaroo 



H=^ + Vo[S^- u(x)] + s„;3(x)V„^[x - u(x)], (4.1) 
2m 



= £'o(k,), (4.3) 
while the gradient correction becomes 

- mSa/3(Xc)[(?)a?)/3) - {Va){v/3)]- (4.4) 



The angular brackets in the above expression represent 
the expectation value of the enclosed operators in the 
Bloch state at k = kc, and Va = dH(k)/dka is the ve- 
locity operator. As is well known, {va) = dSo/dka = Va 
is the group velocity of Bl och electrons. As for the last 
term of the Hamiltonian (4.1), we may write, in accor- 
dance with our approximation [Eq. (2.11)], 



(*|s„;3(i)Va/3|*) ~ S„;5(X,)(V^^), 



(4.5) 



af3 — V'q/3[x — u(xc)]. We again write the energy 



where V, 

of the wave packet in the form 



(4.6) 



where, this time, the correction A£, which is to-be called 
the deformation pot enti al for the wave packet ,l3 has two 
contributions [Eqs. (44) and (4^)]. As the latter is pro- 
portional to the local deformation, it may be written in 
the form 



with 



A£ = S^p (Xe)i^Ji'^(ke 



(4.7) 



(4.8) 



We note that this quantity vanishes in the free electron 
limit (Vo — > 0), which is quite satisfying from a physical 
point of view: a wave packet should not feel the effect of 
a deformation of the lattice in the absence of electron- 
lattice coupling. I— I 

Further, from the differential band structured under 
a uniform strain e , 



£'{k')~£o{k) 



(4.9) 



where k' = ( 1 + s*) ^ • k, A£ can be determined by the 
relation 



Dap = D^p+mVaV/B, 



(4.10) 



according to Eqs. (4.S) and (4.£). The same expres- 
sion has been obtained for the deformation potential for 
electron-phonon interaction at long wavelengths, show- 
ing the equivalence of the |la.tter with our wave-packet 
deformation potential (4.7).e3 



C. Equations of motion and lattice tracking 

Before we proceed further, let us first make the dis- 
placements time dependent, {u;(i)}, as we are discussing 
the dynamical aspect. Accordingly, we extend the re- 
sults of the previous subsections with the replacements 
u(x) u(x, t) and Sq/3(x) — > Sq^(x, t), and with the 
Bloch wave basis given by the states {V'k(x — u(xc, t))}. 



As for the Lagrangian, the functional form of the ba- 
sis states imply that the last two terms of its general 
expression ( p. 19 ) become 



du 



dUa 



du\ dUa 



dt/ 



where 



_md£^ 

- n dk 



hk. 



(4.11) 



(4.12) 



in which h has been restored. The reader is warned that 
the displacement field (u or Ua) should not be confused 
with the periodic part of the Bloch state, \u). The quan- 
tity f (k), being the difference between the group momen- 
tum of a Bloch electron and the momentum of a free elec- 
tron, denotes that part of the momentum arising from the 
lattice interaction alone. It has the desired property of 
vanishing in the free electron limit, where the lattice de- 
formation should not be felt. By substituting Eq. (4.11) 
into Eq. (2.19), we obtain 



L = -£ 



■ kc + kc • 



i^)+u.f(kc), (4.13) 



where u denotes du/dt. The last term is new and repre- 
sents the rate of change of the Berry phase due to lattice 
tracking. We shall see below that there is a tendency for 
the lattice to drag the electron with its displacement mo- 
tion, hence the word "tracking" . This term also gives rise 
to the Burgers vector when integrated around a disloca- 
tion, so in a sense it keeps track of the la ttice position. 
Similarly, the Berry curvatures in Eqs. ( 2.20| ) take val- 



-n 



duj df. 
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^ ° dxcp dk, 

df-y dU-y 



dkrn dt 



(4.14) 



and by expressions ( |4.q ) and (4/7) for the equations of 
motion become 



d£o 

dkrn 



+ S/37 



dD 



dkr 



h — 



Udjca 



ikaXn)a-Ul3-^, (4.15) 

dfp I 



(4.16) 



In the above equations, the semiclassical force has two 
contributions: the first term arises from the deformation 
potential, and the second term, which arises from the lat- 
tice tracking term, will be called the tracking force. For 
uniform strains the semiclassical force vanishes, so that 
kc is a good quantum number in such a case. 

The velocity has three contributions in addition to the 
usual term given by the gradient of the band energy. The 
second term corresponds to the group velocity arising 
from the deformation potential, and the third term, as 
was seen earlier in Eq. (3.8) of section [II, is the Hall 
velocity. The last term arises from the lattice tracking 



term of the Lagrangian, and will be similarly called the 
tracking velocity. To understand this term better, let us 
rewrite the term using Eq. (4.12) as 



Ua — mup 



dkcBdk. 



COL 



(4.17) 



In the free electron limit, where electrons and lattice de- 
couple, the adiabatic velocity vanishes identically, for 



dkcadkcf} 



m 



(4.18) 



On the other hand, when the band under consideration is 



full, the second term of (4.17) averaged over the band is 
zero, implying a complete adiabatic following of the lat- 
tice motion, whickxpufirms an earlier result on adiabatic 
particle transport .cll 



D. Dislocations and Berry phase 

Our formalism is also applicable to dislocation strain 
fields, which are well defined except in a region of a few 
atomic spacings around the line of dislocation. Outside 
this region, the displacement field u(x) is a smooth but 
multiple-valued function. The change in the displace- 
ment field along a closed loop around the line of disloca- 
tion, Au = §^ dxadu/dxa = b, is known as the Burgers 
vector; it equals one of the Bravais lattice vectors. On 
account of this multiple- valuedness, a wave packet of in- 
cident wave vector k taken around the line of dislocation 
acquires a Berry phase 

r = j> dTic-(u = jdvL ■ f (kc) « b • f(k), (4.19) 

where we have assumed {u\i du/dh) to be zero@ and the 
corrections arising from the changes in the wave vector 
to be negligi ble. I n other words the action integral of the 
last term of (4.13) around a dislocation gives us a Berry 
phase that is proportional to the Burgers vector. 

Note that this Berry phase is independent of the path 
as long as it encloses the dislocation line. What wejiiave is 
a situation similar to the Aharonov-Bohm effect ,c2l with 
the dislocation playing the role of the solenoid, and the 
Berry curvature J7xx the role of the magnetic field. Just 
as in the case of the solenoidal field, Cl^x ~ everywhere 
except for the core region where it should be taken to be 
singular. 



The Berry phase (4JJ) affects the scattering of elec- 
trons by a dislocation,^!] and our result can be used to 
compute the shift in the scattering fringe pattern due the 
Berry phase. The intensity distribution has interference 
terms of the form cos(6' -I- F) due to each pair of paths 
along the two sides of the dislocation, where 9 is given 
by the path length difference per wavelength of the in- 
cident beam. The Berry phase and hence the shift are 



maximal when the incident wave vector k is such that 
f (k) = mv — Ti'k is parallel to the Burgers vector b. To fix 
ideas, let us make the assumption that f || k. For an edge 
dislocation, where b lies in the normal plane of the dis- 
location axis, we expect that maximal effect of the Berry 
phase is seen when the direction of the incident beam is 
perpendicular to the axis and coincide with the direction 
of b. For a screw dislocation, the maximal Berry phase 
occurs when k is parallel to the axis along which b lies. 
However, since the beam must pass the dislocation in 
order to produce the interference pattern, the maximal 
Berry phase effect should actually occur when k is along 
some finite angle away from the axis. 

The Berry phase can also affect the electron diffraction 
pattern of a deformed crystal with or without a disloca- 
tion. When an electron beam of wave vector k is sent 
through a crystal, it propagates as Bloch waves in differ- 
ent bands of the crystal, all with the same Bloch wave 
vector which is equal to k modulo a reciprocal lattice 
vector. Because of the energy differences of the bands, 
the Bloch waves grow out of phase from one another as 
they propagate. Therefore, when they exit the crystal, 
they recombine to produce not only the incident beam 
but also the diffracted beams. For a deformed crystal, 
the phase change of a Blo ch w ave is given by the time 
integral of the Lagrangian ( 4.13| ). The contribution from 
the Xc • kc term may be dropped, because it is indepen- 
dent of the band index. For weak strains, we may also 
negle ct th e deformation potential and the third term in 
Eq. (4.13). The dominant correction to the dynamical 
phase arises from the last term: Au • f(k), where Au 
is the change in the displacement field over the path of 
propagation of the Bloch waves. To show the soundness 
of these general ideas, we have calculated the diffraction 
pattern from a thin slab of crystal containing a screw dis- 
location, correctly rseproducing earlier experimental and 
theoretical results. EJ 



V. SUMMARY 



We provide a unified framework for wave-packet dy- 
namics of electrons in slowly perturbed crystals, which 
consists of constructing a wave packet using the Bloch 
states belonging to a single band of the local periodic 
Hamiltonian, obtained by replacing the perturbations by 
their value at the center of the wave packet, and deriving 
its dynamics in a general form, based on the time depen- 
dent variational principle. We derive the wave-packet 
energy up to first order in the gradient of the pertur- 
bations and all kinds of Berry-phase corrections to the 
semiclassical dynamics and the quantization rule. Also, 
we give a discussion of a formal quantization procedure 
through the Hamiltonian formalism with the semiclassi- 
cal dynamics as the starting point. 

We illustrate our framework with two cases of pertur- 
bations. For electromagnetic fields, previous results of or- 



bital magnetization and anomalous velocity are obtained 
purely from a single-band point of view. For deforma- 
tions in crystals, we obtain a Berry-phase term in the 
Lagrangian due to lattice tracking, which gives to the 
new terms of tracking velocity and force in the equa- 
tions of motion of the wave packet. For multiple- valued 
displacement fields in the presence of dislocations, this 
term manifests as a Berry phase, which we show to be 
proportional to the Burgers vector around each disloca- 
tion. Also, we relate the deformation correction to the 
wave packet energy to the shift in band energy under 
uniform strain, which turns out to be the same as the 
deformation potential for electron-phonon interaction at 
long wavelengths. 

The combined effects of electromagnetic and deforma- 
tional perturbations are yet to be studie d. Given that 
the perturbed Hamiltonian is of the form , the equa- 
tions of motion for the wave-packet dynamics can be com- 
pletely determined in terms of the properties of the un- 
perturbed crystal and the fields of perturbations. Such a 
theory will provide a basis for electron transport in de- 
formed crystals and should be pursued in the near future. 
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which are easily verified. 

From Eqs. ( ]AlD , ( |A2[ ), and the completeness relation 



J2 J d^l lV'«q) (V-nql = i, 



(A3) 



we can show that 



(9xc 



■0nq' 



dur, 



{'^P„q\x\^pncl') + 



<5(q-q'). (A4) 



If we assume that the wave packet is constructed from 
the Bloch states of nth band, we have 



(9xc 



(A5) 



q=qc 



as only the second term of Eq. (A2) contribute to the 
expectation value. By the same token, from Eqs. (A5) 
and (A4), we obtain 



axc 



(X - Xc) 



(A6) 



{£c7l 



dUn 



(9q 



q=qc 



Hence the expectation value (2.1S) of Eq. (2.16) is just 
half of the sum of ( A6) and its complex conjugate. 



APPENDIX A: 
THE GRADIENT CORRECTION 

In the derivation below, the band indices for the Bloch 
states and the energy bands have been restored. 

We begin with the matrix elements of the position 



operator,! 

(^/'nqlxlV'n'q') 

and the identity 



. dUn> 

9q 



%'-q), 
(Al) 
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APPENDIX B: 
THE DEFORMED CRYSTAL POTENTIAL 

Given the atomic displacements {u;}, which could be 
due to either strain or rotations, let us denote the de- 
formed crystal potential regarded as a function of the 
coordinates of the electron and the equilibrium atomic 
positions by y(x; {R; -I- u;}). If the displacements are 
small, we can expand the potential in powers of {u;} and 
do away with the perturbation theory. When they are 
large, we have to adopt a different technique. For this 
purpose we introduce a smooth displacement field u(x) 
such that u(Ri -f u;) = uj, and write the potential as 



y[x-u(x);{Rj + uj-u(x)}], 



(Bl) 



where we have used the invariance property of a poten- 
tial under a simultaneous translation of electronic and 
atomic positions. We may now expand the potential in 
powers of u; — u(x) as 



yo[x - u(x)] + ^[u, - u(x)]„V^„ 



(B2) 



where Vq{x — u) = V[x. — u; {R;}] is the potential used 
in the deformable ion model in which the approximated 
potential at x equals the undeformed crystal potential at 
the undeformed electron coordinate x — u(x), and 



Via 



c)y[x-u(x);{R;}] 

dRla 



y/„[x-u(x) -R,]. (B3) 



The equality in the last equation is meant to suggest that 
the center of fall off of Via is at the zero of the expression 
X — u(x) — R/, viz., X = U/ + R;. 



The expansion (B2) is meaningful only if Via decreases 
sufficiently rapidly with increasing |x — u(x) — R/|. This 
condition holds for metals because of Coulorab. screen- 
ing; it has been argued that Via is short rangedES also for 
non-polar semiconductors and insulators. 

Under such a condition, we may write 

Ula = 1ta(x + R; + U; - x) 

W Ua(x) -I- (R; + U; - x)/3Sa^(x) 

W -Ua(x) + [R; + u(x) - x]/3Sa/3(x), (B4) 



whence the summation in Eq. (B2) can be put into the 
form SapVafi with 



Va/3[x - U(x)] = ^[R; + U(x) - x]pVia. 
I 



(B5) 
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